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Quantum  and  Classical  Studies  of  the  0(3P)  +  H2(v=0-3,j=0)  -»  OH  +  H  Reaction  Using 

Benchmark  Potential  Surfaces 

M.  Braunstein*  and  S.  Adler-Golden 
Spectral  Sciences,  Incorporated,  Burlington,  Massachusetts  01 80S 
(♦Electronic  mail:  matt@spectral.com) 

B.  Maiti  and  G.  C.  Schatz 

Department  of  Chemistry,  Northwestern  University,  Evanston,  Illinois  60208-3113 

Abstract 

We  present  results  of  time  dependent  quantum  mechanics  (TDQM)  and  quasi-classical 
trajectory  (QCT)  studies  of  the  excitation  function  for  0(3P)  +  H2(v=0-3,j=0)  — >  OH  +  H  from 
threshold  to  30  kcal/mol  collision  energy  using  benchmark  potential  energy  surfaces  [Rogers  et 
al.  J.  Phys.  Chem.  A  104,  2308  (2000)].  For  H2(v=0)  there  is  excellent  agreement  between 
quantum  and  classical  results.  The  TDQM  results  show  that  the  reactive  threshold  drops  from  1 0 
kcal/mol  for  v=0  to  6  for  v=l,  5  for  v=2  and  4  for  v=3,  suggesting  a  much  slower  increase  in  rate 
constant  with  vibrational  excitation  above  v=T  than  below.  For  H2(v>0),  the  classical  results  are 
larger  than  the  quantum  results  by  a  factor  -2  near  threshold,  but  the  agreement  monotonically 
improves  until  they  are  within  -10%  near  30  kcal/mol  collision  energy.  We  believe  these 
differences  arise  from  stronger  vibrational  adiabaticity  in  the  quantum  dynamics,  an  effect 
examined  before  for  this  system  at  lower  energies.  We  have  also  computed  QCT  OH(v,’j’)  state- 
resolved  cross  sections  and  angular  distributions.  The  QCT  state-resolved  OH(v’)  cross  sections 
peak  at  the  same  vibrational  quantum  number  as  the  H2  reagent.  The  OH  rotational  distributions 
are  also  quite  ‘hot’  and  tend  to  cluster  around  high  rotational  quantum  numbers.  However,  the 
dynamics  seem  to  dictate  a  cutoff  in  the  energy  going  into  OH  rotation.  The  state-resolved  OH 
distributions  were  fit  to  probability  functions  based  on  conventional  information  theory  extended 
to  include  an  ‘energy  gap’  law  for  product  vibrations  and  angular  momentum  constraints  for 
product  rotations. 


1 

Approved  for  public  release;  distribution  unlimited. 


DRAFT  9/16/03 


Quantum  and  Classical  Studies  of  the  0(3P)  +  H2(v-0-3,j-0)  — >  OH  +  H  Reaction  Using 

Benchmark  Potential  Surfaces 

M.  Braunstein*  and  S.  Adler-Golden 
Spectral  Sciences,  Incorporated,  Burlington,  Massachusetts  01803 
(*Electronic  mail:  matt@spectral.com) 

B.  Maiti  and  G.  C.  Schatz 

Department  of  Chemistry,  Northwestern  University,  Evanston,  Illinois  60208-3113 

Abstract 

We  present  results  of  time  dependent  quantum  mechanics  (TDQM)  and  quasi-classical 
trajectory  (QCT)  studies  of  the  excitation  function  for  0(3P)  +  H2(v=0-3,j=0)  — »  OH  +  H  from 
threshold  to  30  kcal/mol  collision  energy  using  benchmark  potential  energy  surfaces  [Rogers  et 
al.  J.  Phys.  Chem.  A  104,  2308  (2000)].  For  H2(v=0)  there  is  excellent  agreement  between 
quantum  and  classical  results.  The  TDQM  results  show  that  the  reactive  threshold  drops  from  10 
kcal/mol  for  v=0  to  6  for  v=l,  5  for  v=2  and  4  for  v=3,  suggesting  a  much  slower  increase  in  rate 
constant  with  vibrational  excitation  above  v=l  than  below.  For  H2(v>0),  the  classical  results  are 
larger  than  the  quantum  results  by  a  factor  -2  near  threshold,  but  the  agreement  monotonically 
improves  until  they  are  within  -10%  near  30  kcal/mol  collision  energy.  We  believe  these 
differences  arise  from  stronger  vibrational  adiabaticity  in  the  quantum  dynamics,  an  effect 
examined  before  for  this  system  at  lower  energies.  We  have  also  computed  QCT  OH(v,’j’)  state- 
resolved  cross  sections  and  angular  distributions.  The  QCT  state-resolved  OH(v’)  cross  sections 
peak  at  the  same  vibrational  quantum  number  as  the  H2  reagent.  The  OH  rotational  distributions 
are  also  quite  ‘hot’  and  tend  to  cluster  around  high  rotational  quantum  numbers.  However,  the 
dynamics  seem  to  dictate  a  cutoff  in  the  energy  going  into  OH  rotation.  The  state-resolved  OH 
distributions  were  fit  to  probability  functions  based  on  conventional  information  theory  extended 
to  include  an  ‘energy  gap’  law  for  product  vibrations  and  angular  momentum  constraints  for 
product  rotations. 
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I.  INTRODUCTION 

Recently,  the  energy-dependent  cross  section  (excitation  function)  for  the  reaction  0(  P)  + 
H2(v=0)  -»  OH  +  H  was  measured  by  molecular  beam  methods  and  was  studied  with  high-level 
quantum  dynamics1  using  ‘benchmark’  potential  surfaces 2  Agreement  was  excellent  over  a  large 
energy  range.  The  above  reaction  is  of  fundamental  interest:  it  is  one  of  the  simplest  hydrogen 
abstraction  reactions,  it  is  important  in  many  combustion  processes,  and  it  is  small  enough  for 
quantitative  theoretical  treatments.  It  has  been  studied  theoretically  with  quantum  '  and 
classical10'12  methods  on  a  number  of  potential  surfaces  and  studied  experimentally  through 
measurement  of  the  rate  constant.13  The  high  barrier  to  reaction  (~13  kcal/mol),  however, 
requires  high  reagent  velocities,  and  this  has  made  the  measurement  of  the  excitation  function 
and  state-to-state  processes  difficult.  The  increasing  capabilities  of  hyper-velocity  O-atom 
sources,  as  were  used  to  study  the  excitation  function,1  should  enable  the  accurate  measurement 
of  state-to-state  quantities  for  this  ‘prototype’  reaction  in  the  near  future. 

In  this  paper  we  present  results  of  time  dependent  quantum  mechanics  (TDQM)  and 
quasi-classical  trajectory  (QCT)  studies  of  the  excitation  function  for  0(3P)  +  H2(v,j=0)  -»  OH  + 
H,  where  the  H2  vibrational  quantum  number  v  varies  from  0-3,  from  threshold  to  30  kcal/mol 
collision  energy,  using  the  same  benchmark  potential  energy  surfaces.  Because  the  vibrational 
spacings  are  large,  quantum  effects  should  be  especially  important.  Comparison  of  the  results 
will  provide  insight  into  the  validity  of  classical  approaches  for  such  a  ‘worst  case.’  In  addition, 
the  present  results  will  enable  us  to  determine  the  change  in  reactivity  with  increasing  excitation 
of  the  reactants,  which  is  a  subject  of  interest  to  fundamental  studies  of  reaction  dynamics.14  In 
past  studies  of  O  +  H2(v=l),  it  was  found  that  the  rate  constant  is  three  orders  of  magnitude 
larger  than  for  O  +  H2(v=0)  at  room  temperature.15  This  implies  a  reduction  in  threshold  energy 
from  the  ~10  kcal/mol  value  observed  for  v=0  in  the  molecular  beam  measurements  to  about  6 
kcal/mol  for  v=l.  If  this  trend  were  to  continue,  then  the  threshold  could  drop  to  zero  for  v=3. 
Large  drops  in  threshold  energy  have  been  found  before,  such  as  for  the  H  +  H20  reaction.14 
However,  the  drop  in  threshold  energy  with  increasing  vibrational  excitation  is  very  sensitive  to 
the  shape  of  the  potential  surface,  so  it  is  not  clear  what  behavior  is  to  be  expected  for  O  +  H2. 
Since  high  quality  surfaces  are  available  for  O  +  H2,  we  can  use  this  reaction  as  a  benchmark  for 

studying  the  drop  in  threshold  with  increasing  vibrational  excitation. 
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An  additional  component  of  the  present  paper  is  that  we  have  computed  QCT  OH(v’,j’) 
state-resolved  cross  sections  and  angular  distributions.  We  show  that  these  OH(v’j’) 

16  17 

distributions  can  be  fit  to  probability  functions  based  on  conventional  information  theory  ’ 
extended  to  include  an  ‘energy  gap’  law'8  for  OH  vibrations  and  angular  momentum 
constraints19,20  for  rotations.  We  compare  all  our  results  to  experimental  data  where  available  and 
to  previous  theoretical  calculations.  These  distributions  will  likely  be  measured  in  future  reaction 
dynamics  experiments,  so  our  results  will  be  useful  for  interpreting  these  measurements. 

The  O  +  H2  (v>0)  reaction  is  also  of  interest  at  high  velocities  in  order  to  characterize  the 
ambient  environment  of  spacecraft  under  low  earth  orbit  (LEO)  conditions.  Spacecraft  in  low 
earth  orbit  continually  emit  gases  due  to  outgassing,  venting,  and  engine  firings.  These  gases  can 
then  interact  with  the  atmosphere,  mostly  oxygen  atoms,  at  high  relative  velocities  usually  near 
~8  km/s.  The  resulting  highly  excited  products,  OH(v’,j’)  in  the  present  case,  can  produce 
radiative  emissions  near  the  vehicle  that  provide  in-situ  monitoring  of  the  spacecraft 
environment. 

The  paper  proceeds  as  follows.  In  Section  II  we  describe  our  computational  methods.  In 
Section  III  we  present  the  detailed  results,  first  focusing  on  the  excitation  functions,  then  the 
state-resolved  product  distributions,  and  finally  the  fits  of  the  distributions.  Section  TV  presents  a 
summary. 

n.  COMPUTATIONAL  METHODS 

Scattering  calculations  were  done  on  the  3 A’  and  3A”  surfaces  of  Rogers  et  al.2  These 
surfaces  have  been  used  previously  to  study  the  O  +  H2(v=0,j=0)  excitation  function  with  good 
agreement  with  measurements.1  The  surfaces  were  generated  from  ab-initio  configuration 
interaction  calculations  with  a  very  large  basis  set,  and  they  are  likely  the  most  accurate  surfaces 
to  date.  Coupling  between  the  two  surfaces  has  been  ignored,  and  the  23A”  state  which  leads  to 
electronically  excited  products  has  not  been  included.  From  previous  studies,  these 
approximations  should  be  of  small  importance.5  Likewise,  any  coupling  to  singlet  surfaces  is 
ignored. 

Results  were  generated  using  3D  time-dependent  quantum  mechanics  (TDQM)  and  quasi- 
classical  trajectory  (QCT)  methods.  The  TDQM  wave  packet  formalism  has  been  described 
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previously,22  including  details  of  the  spatial  grid  and  time  propagation.  To  reduce  the 
computational  burden,  it  has  been  implemented  with  the  centrifugal  sudden  approximation 
(CSA),  which  should  be  a  small  approximation  based  on  previous  studies  of  similar  systems.  For 
the  QCT  calculations,  we  use  standard  Monte  Carlo  techniques.23  The  required  derivatives  of  the 
potentials  were  computed  with  a  fifth  order  finite  difference  method.  At  least  25,000  trajectories 
per  collision  energy  per  electronic  state  per  initial  H2  vibrational  state  were  generated.  For 
H2(v=0),  a  fixed  time  step  of  0.5  a.u.  and  a  maximum  impact  parameter  of  3.0  were  used.  For 
H2(v=1,2)  the  same  time  step  was  used,  but  the  maximum  impact  parameter  was  5.0  a.u.  For 
H2(v=3)  a  time  step  of  0.25  a.u.  and  maximum  impact  parameter  of  7.0  a.u.  were  used. 

III.  RESULTS 

A.  Total  cross  sections 

The  total  cross  sections  for  the  reaction  0(3P)  +  H2(vj=0)  OH  +  H  are  shown  in  Figure  1. 
The  present  results  are  obtained  by  forming  the  statistical  average  of  the  cross  sections  for  the 
two  separate  electronic  states  considered,  (1/3)  (3A’  +  3A”).  For  H2(v=0)  as  seen  in  a  previous 
study,1  there  is  excellent  agreement  between  the  TDQM,  QCT,  and  experimental  results,  except 
at  collision  energies  near  threshold.  Although  the  QM  results  of  Ref.  3  use  the  same  set  of 

potential  energy  surfaces  as  the  present  study,  they  begin  to  diverge  away  from  the  present 

♦ 

results  at  collision  energies  of  15  kcal/mol.  They  are  about  twice  those  of  the  present  results  at 
~23  kcal/mol.  It  is  possible  that  these  differences  come  from  the  ‘J-shifting’  approximation  used 
in  Ref.  3  to  obtain  the  Jtot  >  0  cross  sections.  The  QCT  results  of  Ref.  12  agree  very  well  with  the 
present  QCT  results  at  energies  up  to  ~25  kcal/mol  collision  energy.  For  the  H2(v>0)  cross 
sections  we  shall  see  that  close  agreement  up  to  a  certain  collision  energy  is  also  the  case,  and 
this  may  indicate  that  the  parts  of  the  potential  energy  surface  of  Ref.  12  accessed  at  lower 
energies  may  have  the  same  accuracy  as  the  ones  used  here.  The  QCT  results  of  Ref.  1 1  are 
consistently  too  large,  and  this  is  probably  due  to  inaccuracies  of  the  potential  energy  surface 
used  in  this  study  near  the  barrier  region. 

Examining  the  TDQM  cross  sections  we  see  that  the  reactive  threshold  energy  drops  from 
~10  kcal/mol  for  v=0  to  6  kcal/mol  for  v=l,  5  for  v=2  and  4  for  v=3.  The  v=0  to  v=l  drop  is 
consistent  with  what  would  be  expected  based  on  the  measured  ratio  of  v=l/v=0  rate  constants,15 
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as  mentioned  in  the  Introduction.  In  addition,  this  change  is  consistent  with  expectations 
discussed  for  other  ab  initio  surfaces8  that  there  is  a  change  in  the  location  of  the  reactive 
bottleneck  from  close  to  the  saddle  point  for  v=0  to  the  reagent  side  of  the  saddle  point  for  v=l. 
For  v=2  and  3  we  see  that  the  drop  in  threshold  is  much  smaller  than  that  from  v=0  to  v=l, 
indicating  only  a  modest  shift  in  bottleneck  location  with  reagent  vibrational  excitation.  In 
addition,  the  threshold  remains  finite  for  v=3,  even  though  there  is  plenty  of  energy  to  surmount 
the  barrier. 


O  +  H2(0,0)  ->  H  +  OH  O  +  H2(1 ,0)  H.  H  +  OH 


FIG.  1.  Total  cross  sections  for  0(3P)  +  H2(vj=0)  — »  OH  +  H.  Crosses  are  the  experimental 
results  of  Ref.  1,  solid  lines  are  the  present  TDQM  results,  circles  are  the  present  QCT  results, 
triangles  are  the  quantum  scattering  results  of  Ref.  3,  squares  are  the  QCT  results  of  Ref.  12,  and 
diamonds  are  the  QCT  results  of  Ref.  11. 

For  v>0,  present  quantum  and  classical  results  are  now  only  within  a  factor  of  ~2  near 
threshold  (much  greater  very  near  threshold)  with  the  QCT  results  always  too  large.  As  the  H2 
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vibration  increases  the  agreement  becomes  worse.  The  agreement  monotonically  improves  with 
collision  energy  until  the  QCT  and  TDQM  cross  sections  are  within  ~10%  of  each  near  30 
kcal/mol  collision  energy. 


Figure  2  shows  the  present  TDQM  and  QCT  results  for  0(JP)  +  H2(v,j=0)  ->  OH  +  H,  but 
broken  down  into  contributions  from  each  surface.  As  seen  for  the  total  cross  sections,  the 
agreement  is  very  good  for  H2(v=0),  but  the  QCT  results  are  too  high  for  H2(v>0)  with  closer 
agreement  as  the  collision  energy  increases. 


o  +  H2(0,0)  -►  h  +  OH 


O  +  H2(2,0)  H  +  OH 


O  +  H2(1,0)-»H  +  OH 


FIG.  2.  Cross  sections  from  each  electronic  surface  for  0(3P)  +  H2(v,j=0)  — »  OH  +  H.  Solid  lines 
are  for  the  TDQM  results  using  the  3A’  surface,  broken  lines  are  for  the  TDQM  results  using  the 
3 A”  surface,  circles  are  for  the  QCT  results  using  the  3 A’  surface,  and  squares  are  for  the  QCT 
results  using  the  3A”  surface. 


We  believe  the  differences  between  the  present  TDQM  and  QCT  cross  sections  arise  from 

differences  in  vibrational  adiabaticity  between  quantum  and  classical  mechanics  (with  the 

quantum  results  being  much  more  adiabatic),  an  effect  that  has  been  examined  before  for  this 

system4,8  (and  others24)  at  energies  near  threshold.  As  the  collision  complex  is  formed,  bound 
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motions  perpendicular  to  the  minimum  energy  reaction  path  can  be  defined.  The  quantized 
energy  associated  with  these  bound  motions  added  to  the  minimum  energy  path  create  effective 
energy  thresholds  for  reaction.  If  the  reaction  is  vibrationally  adiabatic,  it  proceeds  along  the 
reaction  path  with  constant  internal  vibrational  energy,  and  it  becomes  sensitive  to  these  effective 
energy  thresholds.  As  has  been  pointed  out,24  due  to  the  uncertainty  principle  these  thresholds 
can  be  important  even  when  the  time  scale  for  the  bound  vibration  is  long  compared  to  the  time 
of  the  reaction.  The  neglect  of  the  quantized  bound  motion  along  the  reaction  path  in  the  QCT 
calculation  can  therefore  make  for  smaller  barrier  heights  than  a  quantum  approach.  This  could 
result  in  higher  QCT  cross  sections,  which  is  presumably  reflected  in  the  present  results.  As  the 
collision  energy  becomes  much  larger  than  the  barrier  heights,  these  effects  should  diminish, 
which  is  what  we  observe.  For  H2(v=0),  the  barrier  location  is  close  to  the  saddle  point  in  the 
potential  energy  surface.  But  for  H2(v>0),  the  effective  barriers  form  in  the  entrance  and  exit 

Q 

channels,  with  the  entrance  channel  bottleneck  expected  to  be  more  significant.  This  may  be 
important  in  understanding  the  difference  between  classical  and  quantum  results  with  increasing 
vibrational  quantum  number.  The  QCT  calculations  also  neglect  other  potentially  significant 
quantum  effects  such  as  tunneling, 

important  close  to  threshold.  0  +  h2(v,o)  ->h  +  oh 


A  number  of  procedures  have  been 
formulated  to  correct  for  vibrational 
adiabaticity  in  QCT  calculations.8,24  A 
simple  approach  is  to  adjust  the  initial 
vibrational  quantum  number  of  the  H2 
to  better  reflect  the  difference  between 
the  collision  energy  and  the  true 
effective  barrier  heights.  We  test  this 
idea  in  Figure  3,  where  we  show  the 
3 A”  TDQM  results  for  v=2  and  the 
corresponding  QCT  results  for  v=l, 
1.5,  1.75  and  2.  Results  for  an  initial 
quantum  number  of  1.5  show  very 


—  TDQM3A"v(H2)  =  2. 
O  QCT  3  A"  VfHj)  =  2. 

□  QCT3A"v(H2)=1.75 
A  GCT3A"v(H2)=1.5 
0  QCT3A'f  v(H2)=1  .0 


FIG.  3.  Comparison  of  TDQM  and  QCT  cross  sections 

for  the  3 A”  surface  alone.  Dashed  line  is  for  TDQM 

with  an  H2  quantum  number  of  2,  circles  are  for  QCT 

with  an  H2  quantum  number  of  2,  squares  are  for  QCT 

with  an  H2  quantum  number  of  1.75,  triangles  are  for 

QCT  with  an  Eb  quantum  number  of  1.5,  and  diamonds 

are  for  QCT  with  an  H2  quantum  number  of  1 .0 
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good  agreement  with  the  TDQM  results  for  most  of  the  energy  range.  However,  we  note  that  the 
OH  product  distributions  are  also  affected,  and  these  show  an  energy  cutoff  which  prevents  the 
higher  (v’j’)  states  from  being  produced,  an  artifact  of  this  simple  procedure. 

B.  State-resolved  cross  sections 

Figure  4  shows  QCT  cross  sections  at  18  kcal/mol  collision  energy  for  the  reaction  0(3P)  + 
H2(v=0-3,j=0)  ->  OH(v’)  +  H,  as  a  function  of  initial  H2  vibration  and  final  OH  vibration. 
Results  at  other  collision  energies  are  similar.  The  results  show  OH  vibrational  distributions 
centered  on  the  same  quantum  number  as  the  initial  H2  quantum  number.  This  makes  for  non- 
thermal  product  distributions  in  which  a  great  deal  of  energy  goes  into  vibration.  (Below  we  fit 
these  distributions  using  an  energy  gap  law  that  depends  on  the  difference  between  the  initial  and 
final  vibrational  energies.)  We  also  show  results  of  QCT  calculations  of  Ref.  12  at  20  kcal/mol 
collision  energy  for  H2(v=0-2,j=l)  where  we  have  multiplied  their  reported  probabilities  by  our 
total  cross  sections  in  order  normalize  out  differences  in  the  overall  reaction  cross  section.  The 
difference  between  the  two  QCT 
calculations  for  H2(v=0)  and  H2(v=2)  and 
j=l  are  so  close  as  to  be  indistinguishable 
on  this  scale.  The  two  QCT  results  for 
H2(v=1)  are  almost  as  close.  Although  at 
this  collision  energy  the  H2(v=2)  total  cross 
sections  differ  by  about  a  factor  of  ~2,  the 
vibrational  distributions  for  all  initial  states 
seem  insensitive  to  any  differences  in  the 
potential  surfaces  of  the  two  QCT 
calculations.  We  also  note  that  Broida  and 
Persky10  have  performed  similar  QCT 
calculations  on  LEPS  surfaces  derived  from 
Ref.  1 1  for  O  +  H2(v=0-7)  ->  OH  +  H.  Like 
the  present  results,  they  found  a  significant 
increase  in  reactivity  with  H2  vibrational 
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QCT  18  kcal/mol 


FIG.  4.  QCT  cross  sections  for  0(3P)  +  H2(v=0- 
3,j=0)  -»  OH(v’)  +  H.  Circles  are  for  H2(v=0), 
triangles  are  for  H2(v=l),  squares  are  for  H2(v=2), 
and  diamonds  are  for  H2(v=3).  Solid  lines  are  for 
the  present  QCT  results  at  18  kcal/mol  collision 
energy,  while  dashed  lines  are  renormalized  QCT 
results  of  Ref.  12  at  20  kcal/mol  collision  energy  for 
H2(v=0-2j=l). 
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excitation  and  an  efficient  conversion  of  vibrational  energy  from  reagents  to  products  (see 
especially  Figure  10  of  Ref.  10).  We  have  examined  several  trajectories  from  our  QCT 
calculations  to  understand  the  origin  of  these  behaviors.  Most  of  the  reactive  trajectories  are 
‘direct,’  with  the  oxygen  atom  leaving  with  the  hydrogen  it  first  encounters  and  the  reaction 
taking  place  in  a  vibrational  period  or  less.  A  quantitative  analysis  of  these  trajectories,  however, 
will  be  required  to  understand  the  mechanism  behind  the  observed  dynamical  behavior.  To 
summarize  our  QCT  results,  Table  1  gives  the  vibrationally  resolved  cross  sections  at  collision 
energies  of  10,  18,  and  26  kcal/mol  and  the  total  cross  sections  for  each  initial  H2(v,j=0)  state. 
Statistical  uncertainties  in  the  QCT  calculations  give  roughly  three  significant  figures  for  the 
overall  reactive  cross  section  from  a  given  initial  H2  vibrational  state,  roughly  two  significant 
figures  for  the  vibrationally  resolved  cross  sections,  and  roughly  one  significant  figure  for  the 
rotationally  resolved  cross  sections  presented  below. 
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Table  1.  QCT  vibrationally  resolved  cross  sections  (summed  over  j’)  in  a.u.  for  0(  P)  + 
H2(vj=0)  ->  OH(v’)  +  H.  The  bottom  of  each  column  shows  the  sum  of  the  v’  cross  sections  as 
well  as  the  TDQM  results  summed  over  v\ 

In  Figure  5  we  show  QCT  results  for  the  rotationally  resolved  cross  sections  for  0(3P)  + 
H2(v=0-3 ,j=0)  -»  OH(v’,j’)  +  H  at  18  kcal/mol  collision  energy.  As  the  initial  H2(v)  quantum 
number  increases  and  more  energy  becomes  available  to  products,  the  OH  rotational 
distributions  peak  at  higher  rotational  quantum  number,  and  for  v’=3  they  correspond  to 
Boltzmann  temperatures  of  -10-20,000  K.  For  H2(v=0)  and  OH(v’=0),  these  results  are 
consistent  with  previous  quantum  studies  (see  Table  V  of  Ref.  5)  done  at  lower  energies  with 
different  potential  surfaces.  As  the  final  vibrational  quantum  state  of  the  OH  increases  and  less 
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3 

energy  is  available  to  rotation,  an  energy  cutoff  becomes  apparent.  This  is  most  evident  for  0(  P) 
+  H2(v=3,j=0)  -»  OH(v’=l-3,j’)  +  H,  where  the  highest  10-15  available  rotational  states  have  a 
negligible  probability.  In  the  next  section,  we  fit  these  rotationally  resolved  probability 
distributions  with  a  functional  form  taking  into  account  the  observed  energy  cutoff. 


OH#)  OH/n 


FIG.  5.  QCT  rotationally  resolved  cross  sections  for  0(3P)  +  H2(v=0-3,j=0)  ->  OH(v’,j’)  +  H  at 
18  kcal/mol  collision  energy,  showing  v’=0  (upper  left),  1  (upper  right),  2  (lower  left)  and  3 
(lower  right).  Circles  are  for  H2(v=0),  triangles  are  for  H2(v=l),  squares  are  for  H2(v=2),  and 
diamonds  are  for  H2(v=3). 

Figure  6  shows  the  QCT  angle  resolved  differential  cross  sections  in  the  center  of  mass  frame 
for  0(3P)  +  H2(v=0-3,j=0)  -»  OH  +  H.  These  cross  sections  are  defined  such  that: 

^ r  (J^coH  ’  J  ~  Q;  _ 1 _ 

~da~  max  N(E,v,j  =  0)  4^(sin^)A^’ 
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where  Nr(E,v,j  =  0,0)  is  the  number  of  reactive  collisions  between  9  and  0  +  A0 .  As  the 


H2(v)  quantum  number  increases,  more 
forward  (near  0  degrees)  scattering  is 
observed.  This  has  been  noted  before  and 
ascribed  to  the  increasing  importance  of  the 
higher  impact  parameters  as  the  vibrational 
quantum  number  increases.  Note  that  for 
v=3,  there  is  a  strong  forward  peak,  indicating 
that  some  fraction  of  the  collisions  at  18 
kcal/mol  involve  a  stripping  mechanism. 


Collision  Energy  1 8  (kcal/mol) 


FIG.  6.  QCT  angle  resolved  differential  cross 

C.  Fits  of  the  state-resolved  distributions 

sections,  dcr/dco  (a.u./sr),  in  the  center  of  mass 

To  father  understand  the  reaction  frame  for  0(3P)  +  H2(v=0-3j=0)  ->  OH  +  H. 

dynamics  and  facilitate  the  use  of  the  H  +  OH  Circles  ^  for  H;(v=0)>  m  for 

product  state  distributions  in  kinetics  TT  ,  ~  tT  ,  , 

v  H2(v=1),  squares  are  for  H2(v=2),  and  diamonds 

modeling,  we  fit  the  distributions  to  simple  _  TT  * 

°  are  for  H2(v=3). 

functional  forms  derived  from  information 

theory.16,17  The  distributions  are  given  in  terms  of  final  state  probabilities,  P;  i.e.,  ratios  of  the 
final  state-resolved  cross  sections  to  the  total  reaction  cross  section  for  a  given  collision  energy 
and  initial  H2  (v,j)  state.  In  information  theory,  probabilities  are  expressed  in  terms  of  a  prior 
(statistical)  distribution  Po  and  a  surprisal  function,  S: 


P  =  P0exp(-S).  (2) 

Simple  formulas  for  Po  are  given  in  the  literature.16,17  For  vibrationally  resolved  distributions,  a 
simple  approximation  gives  Po  °c  (E-Ev>)3/2,  where  E  is  the  energy  available  to  products  and  EV’  is 
the  product  internal  vibrational  energy.  S  is  typically  parameterized  with  a  short  (e.g.,  linear) 
expansion  in  dynamical  variables  such  as  the  vibrational  or  rotational  quantum  number, 
vibrational  or  rotational  energy,  or  energy  between  initial  and  final  states. 
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C.l.  OH  vibrational  state  distributions 

Figure  7  shows  surprisal  plots  (lines  with  symbols)  for  the  OH(v’)  state  distributions  at 
various  collision  energies  with  H2(v=T,3j=0).  The  surprisal  is  defined  -ln(P/Po),  and  it  gives  a 
measure  of  deviation  away  from  purely  statistical  product  distributions.  A  negative  surprisal,  for 
example,  implies  more  energy  in  the  products  than  statistical  considerations  give  alone.  The  “V” 
shapes  of  the  plots  result  from  the  tendency  of  the  reaction  to  conserve  vibrational  energy;  that 
is,  for  the  reaction  probability  to  peak  where  the  OH  vibrational  energy  EV’  equals  the  Ha 
vibrational  energy  Ev.  This  suggests  a  parameterization  of  S  involving  the  vibrational  energy  gap 
|EV-  Ey|.  In  addition,  a  linear  surprisal  parameter  is  needed  to  model  the  low  collision  energy 
cases.  Thus,  the  functional  form  used  to  model  the  vibrational  surprisal  is 

SV’  =  Xv  GV’  +  7.g  |EV’-  Ev|  +  C,  (3) 

where  GV’  is  OH  vibration  energy  relative  to  the  total  product  energy,  (EV’/E).  The  X's  are  fitting 
parameters,  fit  separately  for  each  initial  H2(v)  state,  and  the  constant  C  is  set  by  the  requirement 
that  the  probabilities  over  all  products  for  a  given  initial  condition  sum  to  unity.  All  internal 
energies  are  computed  from  the  spectroscopic  constants  of  Ref.  25  that  are  nearly  identical  to 
those  derived  from  the  potential  surfaces  used. 


FIG.  7.  Vibrational  surprisals  (lines  with  symbols)  and  Eq.  3  fits  (symbols  only)  for  O  +  H2 
(v=l,3)  —>  OH(v’)  +  H  as  a  function  of  the  energy  fraction  in  OH  vibration,  Gv>  =  (EV’/E). 
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The  surprisal  parameters  in  Eq.  (3)  were  determined  from  linear  least-squares  fits  to  the  SV’ 
for  a  given  EC0n  and  v.  When  Xv  and  Xg  are  retrieved  simultaneously,  they  vary  somewhat 
erratically  with  Econ  and  v.  However,  when  an  average  X,g  value  of  0.1 1  per  kcal/mol  is  taken,  the 
retrieved  Vs  behave  smoothly  as  shown  in  Figure  8,  and  the  model  fit  Sv’’s  (as  shown  in  Figure 
7)  still  represent  the  data  well.  The  dependence  of  Xv  on  EC0n  is  approximately  exponential  and 
the  dependence  on  Ev  is  approximately  inverse-square-root.  As  is  commonly  the  case,  the 
vibrational  distributions  become  close  to  statistical  (kv  close  to  zero)  at  high  EC0n-  We  also  note 
that  the  slopes  of  Xy  as  a  function  of  collision  energy  tend  to  zero  with  increasing  H2  vibration. 


O  +  H2(V)  OH  +  H 


FIG.  8.  Dependence  of  X.V)  the  vibrational  surprisal  parameter,  on  the  collision  energy  and  H2(v) 
state. 

A  complication  of  the  QCT  method  is  the  binning  of  product  states  into  discrete  quantum 
state  energy  levels.  If  classical  products  are  near  the  energetic  limit,  the  assigned  quantum  state 
may  be  forbidden  by  energy  conservation.  A  simple  solution  would  be  to  ignore  these 
energetically  forbidden  states.  We  decided  instead  to  include  in  an  approximate  way  some  of  the 
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QCT  information  in  the  fitting  procedure.  This  was  done  by  replacing  a  discrete  quantum  energy 
level  by  a  distribution  of  discrete  energies  that  span  the  classical  energy  bin  for  that  level  up  to 
the  energetic  limit.  (This  yields  OH  fractional  internal  energies  which  can  be  slightly  larger  than 
1.0,  as  seen  in  Figure  7.)  Priors  were  calculated  by  summing  over  this  discrete  distribution.  The 
results  turned  out  to  be  insensitive  to  the  number  of  discrete  points  in  the  distribution  so  long  as 
they  spanned  the  energy  bin.  This  procedure  made  a  modest  change  in  the  fitted  parameters 
compared  to  completely  ignoring  the  energetically  forbidden  quantum  states. 

C.2.  OH  (v’,j’)  state  distributions 

The  information  theory  expression  for  the  OH  vibration-rotation  state  distribution  is16,17 

Pvy  =  Po(v’j’)  (2j’+l)  exp(-S(v\j’)).  (4) 

Here  Po  ~  (Etrans)I/2  =  (E  -  E(v’,j’)) 1/2  is  the  density  of  states  for  free  translation,  where  Etnms  is 
the  product  translational  energy,  expressed  in  terms  of  the  total  product  energy  E  and  the  OH 
internal  energy  E(v’j’).  The  surprisal  S(v’,j’)  can  be  factored  into  vibrational  and  rotational 
parts,  SV’  and  Sy,  which  are  functions  of  vibrational  energy  EV’  and  rotational  energy  Ey, 
respectively.  The  rotational  surprisal,  -ln(Pvy/Po(v\j’))  at  fixed  v’,  may  contain  v’-dependent 
parameters.  Information  theory  suggests  that  the  rotational  surprisal  is  likely  to  be  linear  in  Ey 
and  to  approach  constancy  (statistical  behavior)  at  high  energy. 

Figure  9  shows  plots  of  the  rotational  surprisal  versus  Gr,  the  fractional  available  OH 
rotational  energy,  (Ey/(E-Ev>)),  at  Econ  =  10  kcal/mol  and  26  kcal/mol  for  H2  v=l  and  3.  Like  the 
vibrational  surprisals,  many  of  these  plots  are  also  non-linear.  However,  the  non-linearity,  which 
is  most  pronounced  when  high  rotational  states  are  populated,  is  a  smooth  upward  curvature 
rather  than  a  “V”  shape  characteristic  of  an  energy-gap  law.  This  behavior  is  consistent  with  an 
angular  momentum  constraint,  which  is  not  accounted  for  in  the  conventional  information  theory 
prior  distribution.  We  note  that  for  these  plots  the  normalization  is  such  that  the  Pvy  summed 
over  v’  and  j  ’  is  unity. 


15 

Approved  for  public  release;  distribution  unlimited. 


DRAFT  9/16/03 


O  +  H2(vj=0)  -»  OH(v’,j’)  +  H 
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FIG.  9.  OH  rotational  surprisals  for  H2(v)=l  and  3  using  conventional  information  theory  prior 
distributions.  The  abscissa  is  Gr=  (Ej>/(E-EV’))-  Results  are  given  at  EC0n  =10  and  26  kcal/mol  for 
all  available  OH(v’)  states. 

A  statistical  dynamics  approach  that  incorporates  angular  momentum  constraints  is  the 
PS/OTS  (phase  space/orbiting  transition  state)  model.19  In  this  model,  which  has  been  used  quite 
successfully  to  explain  product  rotational  distributions  from  chemical  reactions  and  unimolecular 
decomposition,  the  precursor  to  the  products  is  an  orbiting  collision  complex  that  adiabatically 
converts  its  rotational  (“orbital”)  energy  into  product  translation.  The  statistical  product 
distribution  is  derived  from  a  direct  count  of  quantum  states  of  the  complex  subject  to 
conservation  of  angular  momentum.  The  product  high-j’  populations  are  constrained  because 
their  high  angular  momentum  must  be  balanced  by  a  correspondingly  high  orbital  angular 
momentum,  which  is  limited  by  the  size  of  the  complex,  defined  by  the  maximum  impact 
parameter  for  the  reverse  reaction.  We  also  note  the  importance  of  the  H  +  OH  product  mass 
combination  in  the  present  case  that  will  limit  the  amount  of  final  orbital  angular  momentum. 
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The  PS/OTS  calculations  are  somewhat  complicated  and  not  very  convenient  for  empirical 
fitting.  As  an  alternative,  we  propose  an  approximate,  analytical  prior  distribution,  Po  ,  which  is 
derived  by  applying  analogous  angular  momentum  and  impact  parameter  restrictions  to  the 
product  translational  density  of  states  (see  the  Appendix).  As  shown  in  Figure  10,  this  modified 
distribution,  which  contains  two  adjustable  parameters,  leads  to  surprisal  plots  that  are  more 
linear  than  those  obtained  with  the  standard  information  theory  Po,  making  it  easier  to 
parameterize  the  QCT  results.  Except  for  a  few  j’  states  near  the  energetic  limit,  the  surprisal 
plots  are  essentially  flat  for  H2(v=3),  indicating  statistical  population  of  the  rotational  states 
within  the  context  of  this  prior  distribution.  We  note  that  priors  incorporating  angular  momentum 
constraints  have  been  used  before,  and  they  have  been  shown  to  greatly  affect  the  internal  state 
distributions.26,27 
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FIG.  10.  As  in  Figure  9,  but  using  the  modified  prior  distributions  with  angular  momentum 
constraints  as  discussed  in  the  text. 
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IV.  SUMMARY 

We  have  presented  results  of  time  dependent  quantum  mechanics  (TDQM)  and  quasi- 
classical  trajectory  (QCT)  studies  of  the  excitation  functions  for  0(3P)  +  H2(v=0-3  j=0)  ->  OH  + 
H,  on  benchmark  potential  energy  surfaces  from  threshold  to  30  kcal/mol.  For  H2(v=0)  there  is 
excellent  agreement  between  quantum  and  classical  results.  The  TDQM  results  show  a  relatively 
large  (4  kcal/mol)  lowering  of  the  threshold  energy  in  going  from  v=0  to  v=l,  and  then  a  much 
smaller  lowering  (about  1  kcal/mol)  in  going  from  v=l  to  v=2  and  from  v=2  to  v=3.  This 
suggests  that  the  effective  bottleneck  for  reaction  moves  up  rapidly  in  energy  in  going  from  v=0 
to  v=l,  but  not  for  higher  v.  In  addition,  this  reaction  provides  an  example  where  even  though  the 
reagent  energy  is  well  above  the  saddle  point  energy,  the  reaction  is  still  activated. 

For  H2(v>0),  the  classical  results  are  larger  than  the  quantum  results  by  a  factor  -2  near 
threshold,  but  the  agreement  monotonically  improves  until  they  are  within  -10%  near  30 
kcal/mol  collision  energy.  We  believe  this  is  due  to  stronger  vibrational  adiabaticity  in  the 
TDQM  results.  QCT  OH(v’,j’)  state-resolved  cross  sections  and  angular  distributions  have  also 
been  computed.  The  QCT  vibrationally  resolved  OH(v’)  cross  sections  peak  at  the  same 
vibrational  quantum  number  as  the  H2  reagent.  The  OH  rotational  distributions  are  also  quite 
‘hot’  and  tend  to  cluster  around  high  rotational  quantum  numbers.  The  state-resolved  OH 
distributions  were  fit  to  probability  functions  based  on  conventional  information  theory  extended 
to  include  an  ‘energy  gap’  law  for  product  vibrations  and  additional  angular  momentum 
constraints  for  product  rotations.  The  fits  were  generally  very  good  for  vibration,  and  the 
modified  surprisals  for  rotation  were  generally  linear  as  a  function  of  collision  energy  for  all 
transitions. 
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Appendix:  Derivation  of  the  modified  prior  distribution  for  OH  v’,j’  states 

We  write  the  prior  distribution  as  the  product  of  the  rotational  and  translational  density  of 
states,  ptrans  Prot-  However,  we  include  only  those  states  satisfying  conservation  of  angular 
momentum.  The  translational  states  included  are  those  whose  reverse  reaction  impact  parameter, 
b,  corresponds  to  an  allowed  range  of  translational  angular  momentum  Lt  =  pbvrei,  p  being  the 
reduced  mass  of  the  products  and  vrei  their  relative  velocity.  At  fixed  vrei  the  density  of  three- 
dimensional  translational  states  is  proportional  to  vrei  and  the  spatial  volume  they  occupy.  In  our 
model  the  volume  is  a  cylinder  with  diameter  br,  the  maximum  impact  parameter  for  the  reverse 
reaction,  oriented  along  the  translation  direction.  Without  any  angular  momentum  restriction,  the 
phase  space  volume  is  proportional  to  the  cross-sectional  area  of  the  cylinder,  %7ibr .  With 
angular  momentum  restricted  to  between  Lmjn  and  Lmax,  the  volume  is  proportional  to  the  annular 

area  %7t(bmax  —  ^min~)  bounded  by  bmax  =  Lmax/pAel  and  bmin  Lmin/pVrei,  where  bmin  —  bmax  —  br. 

Therefore, 

P 0  ~  Prot  Vrel(bmax  ~  bmin  )•  (A  1 ) 

By  angular  momentum  conservation,  the  vector  sum  of  the  collisional  and  reagent  angular 
momenta  equals  the  vector  sum  of  the  product  translational  and  rotational  angular  momenta. 
Therefore,  the  translational  angular  momentum  limits  are  given  by 

Lmax  =  the  smaller  of  Lcon  +  L;  +  Lp  and  pbrvrei  (A2a) 

Lmin  =  the  smaller  of  |  |LC0n  +  Li|  -  Lp  |  and  |  |Lcou  -  Li]  -  Lp  |.  (A2b) 

Here  Lcon  is  the  angular  momentum  associated  with  the  collision  (forward  reaction)  impact 

parameter  and  velocity,  Lj  is  the  angular  momentum  of  the  H2  reagent  and  Lp  is  the  angular 
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momentum  of  the  OH  product.  In  the  present  study  U  =  0.  P0#  may  be  calculated  by  integrating 
over  the  collision  impact  parameter  from  0  to  its  maximum  value.  Alternatively,  a  single  impact 
parameter  close  to  the  RMS  gives  similar  results. 

Specifying  the  magnitudes  of  Lcon,  Lp  and  Lt  uniquely  specifies  the  direction  of  the  OH 
angular  momentum  vector.  Therefore,  only  a  specific  mj  component  of  the  j  manifold  should  be 
counted  in  the  rotational  density  of  states.  As  a  result,  prot  must  contain  the  factor  l/(2j’+l), 
which  cancels  the  2j  ’+1  factor  appearing  in  Equation  (4). 

Equation  (Al)  depends  on  the  collision  energy,  product  translational  energy,  and  OH 
rotational  energy  as  well  as  on  the  forward  and  reverse  reaction  maximum  impact  parameters. 
For  long-range  R  ^  potentials,  the  maximum  impact  parameter  scales  as  Etrans  (see  Ref.  18). 
Combining  constants  and  setting  Li  =  0  leads  to  the  simple  functional  form 

P0#  ~  Etrans1 /6(fmax  -  fmin)/(2j’+l)  (A3) 

where 

fmax  =  the  smaller  of  1  and  (kiEcoii1/3  +  k2Ej1/2)2/Etrans2/3  (A4a) 

fmin  =  the  smaller  of  1  and  (kiEC0n1/3  -  k2Ej1/2)2/Etrans2/3-  (A4b) 

In  Eq.  (A4),  ki  and  k2,  which  are  related  to  the  maximum  impact  parameters,  are  taken  as 
empirical  fitting  parameters.  By  trial  and  error  we  have  found  suitable  values  for  the  H2  +  O 
reaction  to  be  k]  =  0.2  and  k2  =  0.17  kcal/mol'1/6.  Similar  P0#  results  are  obtained  with  forward 
impact  parameter  integration,  in  which  a  weighted  averaging  of  Eq.  (A3)  is  performed  over  the 
k]=0  to  0.3  range. 
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